Role of disorder in half-filled high Landau levels 
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We study the effects of disorder on the quantum Hall stripe phases in half-filled high Landau 
levels using exact numerical diagonalization. We show that, in the presence of weak disorder, a 
compressible, striped charge density wave, becomes the true ground state. The projected electron 
density profile resembles that of a smectic liquid. With increasing disorder strength W, we find 
that there exists a critical value, W c ~ 0.12e 2 /ef, where a transition/crossover to an isotropic phase 
with strong local electron density fluctuations takes place. The many-body density of states are 
qualitatively distinguishable in these two phases and help elucidate the nature of the transition. 
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PACS numbers: 73.20Dx, 73.40.Kp, 73.50Jt 

Recent experiments have discovered a set of new, com- 
pressible states with anisotropic low-temperature mag- 
netotransport properties in high mobility quantum Hall 
systems near half-filling of the high Landau levels (LLs) 
[|l| . Subsequently, it was proposed that these properties, 
emerging as a consequence of strong long-range electronic 
correlations and the quenching of the kinetic energy, are 
manifestations of the electronic liquid-crystal phases with 
partially broken translation and rotational symmetries 
||,||. The anomalous anisotropic transport can be inter- 
preted as arising from a rotational symmetry breaking 
ground states, i.e. a smectic or a stripe ordered charge- 
density wave (CDW) phase @§@jU|. Remarkably, a 
unidirectional CDW state has been predicted to occur 
in half-filled high LLs in the Hartree-Fock (HF) theory 
0. However, it has been known that the smectic liquid 
phase in the HF theory is always unstable to translation- 
symmetry breaking density modulations along the stripes 
|||,||J|J|,P . The latter produces an anisotropic but incom- 
pressible striped crystal that is insulating in both direc- 
tions, inconsistent with experiments. It is thus important 
to consider quantum fluctuations beyond the HF theory, 
which is under active investigation. 

In this paper, we study the effects of disorder on the 
ground state and low energy excitations in half-filled high 
LLs using exact numerical diagonalizations of finite sys- 
tems with up to twelve electrons. Previously, exact di- 
agonalizations were carried out by Rezayi, Haldane, and 
Yang H (hereafter referred to as RHY) for clean sys- 
tems. Although the ground state in this case does not 
directly exhibit stripe order, RHY discovered that the 
static density response function x(l) an d the correla- 
tion function So (q) are both sharply and strongly peaked 
at q* = (q* , 0) , indicative of the tendency towards the 
formation of a unidirectional CDW. They conjectured 
that in the thermodynamic limit, an external modula- 
tion with a period matching q* will result in a quantum 
Hall stripe phase. We will show that the presence of dis- 
order leads to several interesting and important results 



for the ground state and the low energy excitations: (i) A 
very weak random potential associated with disorder in- 
duces a compressible quantum smectic phase as the true 
ground state, (ii) There exists a critical/characteristic 
disorder strength, W c ~ 0A2e 2 /e£ where £ is the mag- 
netic length, which marks a transition/crossover from the 
anisotropic smectic phase at W < W c to an isotropic fluid 
phase at W > W c ||. We illustrate the changes in the 
collective excitation spectrum in terms of the disorder- 
averaged many-body density of states (mDOS) at low 
energies. We show that the latter exhibits an intrigu- 
ing crossover from that resembles a finite density of in- 
gap states in the smectic phase to that of a vanishing 
Coulomb gap behavior in the isotropic phase. 

We consider a two-dimensional system in an a x b rect- 
angular cell with periodic boundary conditions imposed 
in both x and y directions fl(i| . In the presence of a strong 
magnetic field, one can project the Hamiltonian onto the 
up-most, partially-filled, N-th Landau level. The pro- 
jected Hamiltonian in the presence of both Coulomb in- 
teraction and disorder can be written as: 
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where is the guiding center coordinate of the i-th 
electron ||, L^{x) is the Laguerre Polynomial, V{q) = 
2ire 2 /eq is the Coulomb potential, and Vi m p(g) is the im- 
purity potential. We set 1 = 1 and e 2 /el = 1 for con- 
venience. In Eq. ([jj), the wave vector q = 2iz(i x /a,i y /b) 
with i x and i y integers. The disorder potential is gen- 
erated according to the correlation relation in g-space 

< Vi mp (q)Vi mp (— q') >— W 2 5 q - q >, which corresponds to 

< V(r)V(r') >= W 2 S(r — r') in real space, where W is 
the strength of the disorder in units of e 2 /e£. We con- 
sider up to N e — 12 electrons, spanning a Hilbert space 
of size A^basis = 2, 704, 156. Since the disorder potential 
spoils the symmetry classifications of the states used in 
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RHY, we must deal with the entire Hilbert space. We 
obtain the exact low energy eigenstates and eigen wave- 
functions using the Lanczos diagonalization method. 
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FIG. 1. The projected electron density p(r) at four dif- 
ferent disorder strengths for N e = 10, N = 2 and aspect 
ratio A = 0.75. (a) W = 0.015. (b) W = 0.045. (c) 
W = 0.10. (d) W = 0.14. 

The most direct probe of charge density order in the 
ground state is the projected electron density p(r), i.e. 
the equivalent electron density in the lowest LL describ- 
ing the spatial distribution of the guiding center (GC) 
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It is given by 



J'1,3'2 



(2) 



where the summation is over all single particle orbits and 
4>j(r) is the corresponding wavefunction in the lowest LL. 
In Fig. 1, we show how p(r) in the ground state evolves 
with disorder strength W for the number of electrons 
N e = 10 in the N = 2 LL. The aspect ratio of the sys- 
tem is A = b/a = 0.75, which is an optimized geometry 
for observing the instability toward the formation of a 
stripe phase as suggested by RHY. For a very weak dis- 
order strength, W = 0.015, the typical behavior of p(r) 
is shown in Fig. la as a function of the two-dimensional 
coordinates. Two nearly perfect stripes are formed along 
the y-direction with essentially no density modulations 
along the stripes. This is the first direct demonstration 
by exact method of the stripe ordered ground state in 
finite systems in the presence of weak disorder. 

It is important to note that in the clean case studied by 
RHY, the ground state and all exact many-body eigen- 
states are spatially homogeneous due to the translational 
symmetry of the Hamiltonian. Only certain combina- 
tions of the five low- lying, nearly degenerate eigenstates, 
separated by the unique wave vector q*, can give rise to 
a stripe phase that has a substantial overlap with the HF 
unidirectional CDW state. It may be possible to realize 



such a constructed stripe phase in the thermodynamic 
limit in clean systems if an external potential modula- 
tion is present with a wave vector that matches q* [0. 
However, when a relatively weak disorder is turned on, it 
causes a natural mixing of the low-lying eigenstates and 
the near degeneracy is lifted. At the same time, the ran- 
dom disorder potential necessarily introduces an oscilla- 
tory component matching the wave vector q* . The pres- 
ence of the enormous static density response x(q - * q*) 
then gives rise to the unidirectional CDW, or the stripe 
order in the ground state that we observe. 




FIG. 2. The static density-density correlation function 
So(q) at (a) weak disorder W = 0.09 < W c and (b) strong 
disorder W = 0.14 > W c . 

To further support the picture of weak disorder in- 
duced stripe order in the ground state, we calculated the 
projected equal time density-density correlation function 
5o(q) defined as follows: 



So(q) = ^<0|5>W- R ; 



] \0 > 
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We find that for weak disorder of strength W < 0.05, 
S'o(q) is quantitatively indistinguishable from the result 
of RHY for the corresponding clean system, showing 
strong and sharp peaks at q* = (±0.97,0) for 2n/a = 
0.485, with a peak value So(q*) = 3.7 in the case of 
N e = 10. The peak value of 5o(q) increases with the 
number of electrons to 4.9 at N e — 12. It can be seen 
from Fig. la that the width of the stripes is given by 
D s = a/2, which equals 2ir/q* with q x — Air/a, as im- 
plied in this picture. It is interesting to note that the 
edges of the stripes are rather soft, resembling that of 
a single dominant harmonic variation. Next we discuss 
the evolution of the ground state with increasing disor- 
der strength W. Fig. lb shows the real space profile of 
p(r) at W = 0.045. There are still two complete stripes 
with the same stripe width D s . Although density mod- 
ulations along the stripes are still minimal, the shape of 
the stripes exhibits pronounced variations in the stripe 
direction. Thus we conclude that for weak to moderate 
disorder strength, the ground state is in the quantum 
Hall smectic phase 

The behavior of the static density-density correlation 
function Sb(q) m the smectic phase is exemplified in 



2 



Fig. 2a for W = 0.09. Despite the relatively strong dis- 
order potential, the peaks of So(q) are still locked at q*, 
but the peak value has been reduced to 1.96. The en- 
tire structure of S'o(q) remains qualitatively unchanged 
from the clean limit S. The absence of new, emerging 
peaks at other q suggests that the smectic ground state 
is stable against density modulations along the stripes, 
even in the presence of translational-symmetry breaking 
disorder potentials. 

We next demonstrate that further increasing the dis- 
order strength W eventually makes the smectic phase 
unstable. The projected density p(r) is shown in Fig. lc 
for W = 0.1. At this value of W, the smectic phase is 
on the brink of being unstable. p(r) shows very strong 
shape fluctuations and significant density modulations 
along the stripes with detectable inter-stripe correlations. 
At a still larger W = 0.14, the p(r) plotted in Fig. 2d 
shows that the long-range stripe order is lost: the stripes 
become raptured and riddled with defects, very much 
like in the nematic phase. Furthermore, the local ori- 
entation of the short-range stripes show indications that 
it has been rotated away from the y-axis. The corre- 
sponding So(q) is shown in Fig. 2b. It is clear that the 
primary peaks around q* become very much broadened 
with significantly reduced peak values. Concomitantly, 
several competing peaks appear in So(q), marking the 
breakdown of the smectic phase by random disorder and 
the emergence of a nematic/isotropic phase with local 
electron density variations. The positions of these peaks 
are random and dependent on the disorder configura- 
tion. These results suggest that there exists a smectic 
to nematic-like or isotropic phase transition, or a sharp 
crossover that takes place at a critical/characteristic W c 
between W = 0.10 and W = 0.14. We note that it is 
somewhat difficult to distinguish between a nematic-like 
and an isotropic phase with short-range stripe correla- 
tions in our finite size systems. 

In order to pin down W c , 40 different disorder config- 
urations have been studied in detail. We find that the 
disorder strength W c , at which multiple random peaks 
develop in Sq (q) and compete with the primary peaks at 
q*, is always around a narrow region of W c = 0.12. As 
expected, the same value of W c can be extracted directly 
from the spatial distributions of p(r) shown in Fig. 1. 
In Fig. 3a, we plot the disorder averaged peak values of 
5 (q*) for N e = 8 and N e = 10 as a function of W. 
So(q*) stays almost constant for relatively weak disor- 
der (W < 0.05). Then it begins to drop rather quickly 
around W — 0.07, while the strong peaks of Sb(q) remain 
at q* and the emergent peaks at other values of q are ab- 
sent. For even stronger disorder, 5o(q*) becomes compa- 
rable to the background value of S*o(q), which is approxi- 
mately 0.5 for electrons at half-filling, and simultaneously 
random peaks emerge around q*. Such a change occurs 
around W c — 0.12 for N e = 10 as marked in Fig. 3a. 
Due to the finite size effect, the value of W c for N e = 8 



is around 0.1. 
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FIG. 3. (a). The peak vdlme of the So(q) at the primary 
wave vectors q = q* as a function of W with a critical 
value W c — 0.12. (b) The constructed "order parameters" 
of the smectic phase A x 



and A y as a function of W. 



We next construct the "order parameters" of the smec- 
tic phase A x and A y , which essentially describe the over- 
lap of the true ground state with the HF unidirectional 
CDW state fll, 



=< 



(4) 



Here, < ... > denotes the average over 200 and 40 dis- 
order configurations for N e = 8 and N e = 10 respec- 
tively. In Eq. (^), the summation is over all single- 
particle cyclotron orbits along x-direction with a density 
rij x =< 0|a^aj x |0 >, n s = 2N e is the number of orbits 
in each LL whereas jo = (D s /2)/(a/n s ) = n s /4 is the 
half stripe width in units of orbit number. A v can be 
similarly defined. 

In Fig. 3b, both and A y are shown as a function of 
the disorder strength W. At W = 0, the ground state is 
in the uniform phase for the translational invariance of 
the system, A x = A y = 0. However, when we switch on a 
weak disorder potential that stabilize the stripes stacked 
in the x-direction, becomes much larger than A y and 
is close to 1. Only for the HF wavefunction, where the 
cyclotron orbits are either fully occupied or completely 
empty, A^ F = 1. The deviation of A^, from 1 at small 
W is due to the smooth variation of < 0|a^"Oj|0 > be- 
tween high and low occupation orbits, i.e. the softness of 
the edges shown in p(r). This remains true for N e = 12. 
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Although larger system sizes are needed to rule out the 
finite size effect, we point out that the absence of the 
sharp edges is an important observation since many of 
the recent theoretical descriptions of the smectic phase 
^HHH s h are a common starting point of coupled onc- 
dimensional edge states along the stripes, which requires 
the presence of sharp edges. 

When W is increased to 0.1, A^ drops quickly to about 
0.3, corresponding to the development of the local density 
modulations along the stripes shown in Fig. lc. Again 
we see that W c ~ 0.12 can be considered as a transi- 
tion point between a smectic and an isotropic phase with 
comparable A x and A y , signaling the loss of long-range 
anisotropic charge ordering in the ground state. In the 
strong W region, both and A y effectively describe the 
strong local charge density fluctuations. 




FIG. 4. The mDOS verses energy E at three different 
values of W for iV e = 10 and N = 2. 



To develop further insights into the physical origin 
and the nature of such a smectic, quantum Hall liquid 
state, we plot, in Fig. 4, the many-body density of states 
(mDOS) for N e — 10 in the second LL. It is clear that for 
weak disorder, represented by W — 0.015, the mDOS ex- 
hibits features of a finite density of in-gap states. Specif- 
ically, there is a finite density of low-lying states sepa- 
rated from all excited states by a finite energy gap. This 
is a direct evidence that the ground state is compress- 
ible. The gap value is E g ~ 0.015 for both N e = 10 and 
N e = 8. All states below the gap (its number equals N e ) 
have a similar So(q) and projected electron density p(r), 
each represents a stripe state in the y-direction shifted 
away by one orbit distance. They correspond to long- 
wavelength phonon excitations. The higher energy states 
above the gap are isotropic, corresponding to interstitial 
and defect-like excitations accommodated by the impu- 
rity potential. The presence of such a finite energy gap in 
weak disorder and the nature of the in-gap states suggest 
that the smectic state obtained here is likely to be the 
stable ground state in the thermodynamic limit. 

Upon increasing W , the isotropic states move toward 
lower energy and the gap is reduced. At W = 0.045, 
the gap region has turned into a minimum in the mDOS. 
However, the mDOS remains finite in the limit E — > 0, al- 



though its value has been reduced as states are transfered 
to higher energies by disorder. The filling of the energy 
gap begins to have an impact on the projected electron 
density in the ground state shown in Fig. lb in the form 
of pronounced shape fluctuations along the stripe direc- 
tions. With the further increase of W, the gap eventually 
collapses. At W — 0.09, the mDOS as a function of E 
exhibits a remarkable behavior: it vanishes linearly in 
the limit E — > 0. This behavior remains qualitatively 
unchanged in the strong W region and is insensitive to 
sample sizes for N e = 8, 10 and 12. Further increase of 
W in the absence of the spectral gap causes strong mix- 
ing of the low-lying and the excited states and eventually 
destroy the smectic phase in favor of an isotropic phase 
with short-range stripe correlations. 

We have also studied the transport property of the 
anisotropic CDW state by calculating the Thouless en- 
ergy AE XX or AEyy, which is the change in the ground 
state energy under the change of the boundary condition 
in the x or y direction from periodic to antiperiodic. We 
find that for weak disorder, e.g. W = 0.015, AE yy is 
about 50 times bigger than AE XX , suggesting |Q a large 
transport anisotropy associated with the quantum Hall 
smectic phase consistent with the experimental findings. 

We would like to acknowledge helpful discussions M. 
Fogler, E. Rezayi and M. P. A. Fisher. D.N.S. would 
like to thank the hospitality of the Aspen Center for 
Physics where this work was initiated. D.N.S. is sup- 
ported by Petroleum Research Fund grant No. ACS- 
PRF # 36965- AC5 and Research Corporation award No. 
CC5643. Z.W. is supported by DOE grant No. DE- 
FG02-99ER45747 and an award from Research Corpora- 
tion. 



M. P. Lilly, K. B. Cooper, J. P. Eisenstein, L. N. Pfeif- 
fer, and K. W. West, Phys. Rev. Lett. 82, 394 (1999); R. 
R. Du, D. C. Tsui. H. L. Stormer, L. N. Pfeiffer, K. W. 
Baldwin, and K. W. West, Solid States Commun. 109, 389 
(1999). 

E. Fradkin and S. A. Kivelson, Phys. Rev. B 59, 8065 
(1999). 

A. H. MacDonald, and M. P. A. Fisher, Phys. Rev. B 61, 
5724 (2000). 

H. A. Fertig, Phys. Rev. Lett. 82, 3693 (1999). 

A. L opatnikova, B. I. H aperin, S. H. Simon, and X.-G. 



[2] 

[3] 

[4] 
[5] 

[6] 

[7] 



[9] 



:ond-mat/0 105079 



Wen, 

D. G. Barci, E. Fradkin, S. A. Kivelson, V. Oganesyan, 
Phys. Rev. B65, 245319 (2002). 

A. A. Koulakov, M. M. Fogler, and B. I. Shklovskii, Phys. 
Rev. Lett. 76, 499 (1996); R. Moessner and J. T. Chalker, 
Phys. Rev. B 54, 5006 (1996). 

E. H. Rezayi, F. D. M. Haldane, and K. Yang, Phys. Rev. 
Lett. 83, 1219 (1999); it ibid 85, 5396 (2000). 

Note that the melting temperature of the stripes in 
the Hartree-Fock theory vanishes at a critical disorder 



4 



strength. T.D. Stanescu, I. Martin, and P. Phillips, Phys. 

Rev. Lett. 84, 1288 (2000). 
[10] All results remain essentially unchanged under antiperiodic 

boundary conditions. 
[11] For reviews, see e.g., The Quantum Hall Effect, edited by 

R. E. Prange and S. M. Girvin (Springer, New York, 1990), 

2nd ed. 

[12] D. G. Barci, E. Fradkin, Phys. Rev. B65 245320 (2002). 



[13] M. M. Fogler, fcond-mat/01 11001 . 

[14] For a noninteracting system, Thouless energy AE XX or 
AE yy can be related to the conductance quantitatively. 
For an interacting system, a xx should still be a monotonic 
function of AE XX as it measures the general extensiveness 
of the wavefunction by its sensitivity to the change in the 
boundary condition. 



5 



